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o 

J>~j' Abstract. For the isotropic Lame system, we prove in dimensions three or larger that 



both Lame coefficients are uniquely recovered from partial Cauchy data on an arbitrary 
open subset of the boundary provided that the coefficient [i is a priori close to a constant. 



In a bounded domain Q C R , iV > 3 with smooth boundary we consider the isotropic 



& ' Lame system: 



(0.1) 7^.( C ^7^)=° infi,l<z<iV 



],k,l=l J 

and 

(0.2) u\ dn = f, 

^ ■ 

■ where 
CO ' 

Cijki = \{x)5 ij 5 k i + ii(x)(S ik Sji + S u 5 jk ), 1 < i, j, k,l <N 

where the Kronecker delta is denoted by 8ij. The functions A and /i are called the Lame 
q \ coefficients, u(x) = (ux(x), . . . ,un(x)) is the displacement. Assume that 

(0.3) n(x) > onH, (A + /i)(x)>0 on H. 



We set 



£ : - 1 2. ^-«^-' ^« ^ 

\j,k,l=l ' j,k,l=l 

where v — (yi, . . . , vn) is the outward unit normal vector to dfl. The map Ax^u is called 
the Dirichlet-to- Neumann map that maps the displacement at the boundary to the traction. 
Denote 

The partial Cauchy data is defined by 

' — ' 3 

C\,p = {(u, A Xlfl u)\f, Cx^(x,D)u = infl, «|an = /, supp / c T, / e # 5 (dQ)}. 

Here T is an arbitrarily fixed open subset of dQ. We set r = dQ \ T. 

In this paper, we consider the following inverse problem: Suppose that the partial Cauchy 
data C\ jfl are given. Can we determine the Lame coefficients A and \x ? 
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This inverse problem has been studied since the 90 's. A linearized version of this inverse 
problem for full data was studied in [Bj. In two dimensions, Akamatsu, Nakamura and 
Steinberg P proved that for the case of full Cauchy data (T = dfl) one can recover the Taylor 
series of the Lame parameters on the boundary provided that the Lame coefficients are C°° 
functions. This boundary determination result was extended in [13] to higher dimensions. 
In [11] Nakamura and Uhlmann, for the case of full Cauchy data, established that in the two 
dimensions the Lame coefficients are uniquely determined, assuming that they are sufficiently 
close to a pair of positive constants. Recently Imanuvilov and Yamamoto in [9] proved for 
the two dimensional case that the Lame coefficient A can be recovered from partial Cauchy 
data if the coefficient \i is some positive constant. For the three dimensional case Nakamura 
and Uhlmann in [12], [T3] and independently in [5] Eskin and Ralston proved uniqueness 
for both Lame coefficients provided that /i is close to a positive constant. The proofs in 
the above papers rely on construction of complex geometric optics solutions. On the other 
hand, unlike the case of the Schrodinger operator, for partial Cauchy data, the construction 
of such a solutions for the isotropic Lame system seems to be possible only for a dense set 
of Lame coefficients. 

The recovery of Lame coefficients by partial Cauchy data on an arbitrary subboundary is 
useful from the practical point of view, because one can limit input and measurement subsets 
of dQ as much as possible. To the best of our knowledge, there are no results on the unique 
recovery of the Lame coefficients from the partial Cauchy data in the three dimensional case. 
The purpose of this article is to prove such uniqueness in three dimensions. 

Finally we mention that this inverse problem is closely related to the method known as 
Electrical Impedance Tomography (EIT). EIT is used in prospection of oil and minerals and 
in medical imaging in detecting breast cancer, pulmonary edema, etc. For the mathematical 
treatments of this problem, we refer to [2], [3], [I], [7], [8], [10], [15] and the review paper 



Our main result is the following theorem. 

Theorem 0.1. Let fi2 be some positive constants and X±, A 2 G C°°(Q) be some functions 
satisfying (GO)) and X 1 = A 2 on T . If ' C Xl , ill = C A2iM2; then (Ai,/^) = (A 2 ,/i 2 ). 

Proof. The proof consists in showing that from partial Cauchy data one can recover the 
full Cauchy data. First following |13] we obtain that 



1161. 



(0.4) 




Let Ui e H 2 



(Q) be a functions such that 



(0.5) 





where supp/ C T. Since the partial Cauchy data are the same, we obtain 



(0.6) 



Aa 1jW wi = A A2iM2 w 2 on T. 
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This equality combined with ( 10. 4 1) implies that 

(0.7) — ) = (« 2 , — ) onT. 

Since the functions /ij are assumed to be constants, from ( 10.41) we conclude that 
(0.8) Hi = H2 on Q. 

For constant m, we note that 

C\ tfl (x, D)u = fiAu + (/i + A)Vdivit + (divw)VA. 
Applying to equation (10. 5p the operator rot and using the fact that \ij is constant, we obtain 
(0.9) fijArotUj = in Q. 

From ( 10 .4p . ( 10. 7p and equation (10. 5p we conclude 

c)u c)u 
(0.10) {ui,—,d XiXk u x ) = {u2,-^,d XiXk u 2 ) on f , Mi, fc G {1, 2, 3}. 

Hence 

<9rot ?/i cfrot u 2 ~ 

(0.11) (rot Mi, ) = (rotn 2 , ) on T. 

Equality ( 10. lip and the uniqueness of the solution for the Cauchy problem for the Laplace 
equation imply 

(0.12) rot Mi = rotM 2 in Q. 

The Lame operator, with the coefficient « = const, can be written in the form C(x, D)u = 
V((A + 2u)divM) — urot rotw. Then using ( IQ.Sp . ( I0.12p we obtain 

(0.13) V((Ai + 2/ii)divui) = V((A 2 + 2^i)divM 2 ) mil. 

Since (Ai + 2/ii)diviti = (A 2 + 2//i)divw 2 on T, equation ( I0.13P implies 

(0.14) (Ai + 2 / ui)divMi = (A 2 + 2/ii)div-u 2 in Q. 

From (EHD, ( 1TU41) . (10121) and the assumption (Ai - A 2 )|r = we conclude 

dui duo 

(0.15) ^ = ^ onr . 

Therefore if supp / C f in fLl)and / G then 2a = *ia on an. 

Next let / G H^(dil) and the functions G H 2 (Q,) be solutions of the following boundary 
value problem 

(0.16) C Xj , H (x,D) Vj = mQ, v j \ dn = f, j e {1,2}. 

We claim that 



< 017) * a. 



dvi dvx ~ 
-n — = — on 1 . 



4 



O. IMANUVILOV, G. UHLMANN, AND M. YAMAMOTO 



Indeed, let Wj G H 2 (Q) be a solution to the Lame system 
(0.18) C\ jtH (x,D)w j = infi, Vj\ m = g, j G {1,2}, 

3 ' — ' 

where g G H^{dVt) and supp g C T is an arbitrary function. Taking the scalar product of 
equation (10. 16f) with Wj and integrating by parts, we have 

= / {Cx^^x.D^j.w^dx = / (v j ,C Xj , N (x,D)w j )dx + / ((A^tfj, u/j) — (A\^Wj, Vj))d(T 
Jn Jn Jan 

((A Xj ,^Vj,g) - (Axj^wj, f))da = l(A Xjjflj v j: g)da - / (A Xj , H Wj, f)da 
an Jr Jan 



L( A ^,H v i>9)d(T - / (Ax^W!, f)da, 
Jr Jan 



where da denotes the surface measure. 
This integral identity implies 

Aa 1iMi Vi = A X2:fl2 v 2 on f . 
Repeating the arguments fl0.9p - fl0.15p we conclude 

(0-19) -^ = ^ onr . 

Hence, by f l0.15p . (10. 19j) the following full Cauchy data are equal: 

where 

' — ' 3 

= {(u,A x ^u)\ dn ; C x ,n(x,D)u = in fi, u\ dn = f, f e H2(dtt)}. 

Applying the result of [5], [12] and [H], we obtain that Ai = A 2 . □ 
This result immediately implies a local result for \x near constant. 

Corollary 0.1. Let B be a bounded set in C°°(Q), Ai, A 2 £ B, X± — A 2 on T and Xj(x) > 
C > 0,fij(x) > C > on Q. There exist positive e(B) > and positive sufficiently large 
number N such that if^2 k=1 \\ ^f^k\\c N (n) — e (-^) an dCx ltfll = Cx 2 ^ 2 , then (Ai,/xi) = (A 2 ,/i 2 ). 

Proof. Our proof by contradiction. Suppose that the statement of the corollary is false. 
Then there exists a sequence of positive {ej} c ^ =1 such that €j — > and for each ej there exists 
{(Afej, /ifcj)} such that 

2 



(0.20) {h,j}T= l ^B, k = {1, 2}, and £ || V// fe „ 



k=l 



and 

(0-21) C Ali ., Mli .=C A2JjM2i . Vje{l,...,oo}. 

By (10.201) and (10.211) there exist Ai, A 2 G C°°(f2) and positive constants Hk such that 



Applying Theorem I0.1[ we complete the proof of the corollary. □ 
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